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The CoGeNT collaboration has reported evidence of an annual modulation in its first fifteen 
months of data. Here we compare the amplitude and phase of this signal to the modulation ob- 
served by the DAMA collaboration, assuming that both arise due to elastically scattering dark 
matter (DM). We directly map the CoGeNT signal to the DAMA detector without specifying any 
astrophysical parameters and compare this with the signal measured by DAMA. We also compare 
with constraints from CDMS II and XENONIO. We find that DM of mass 5-14 GeV that cou- 
ples equally to protons and neutrons is strongly disfavoured. Isospin-violating DM fares better but 
requires a boosted modulation fraction. 
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I. INTRODUCTION 

There is strong observational evidence for a large abun- 
dance of particle dark matter (DM) in our Universe. 
Detecting this matter through non-gravitational interac- 
tions is a great challenge. A potentially characteristic 
signal arising from particle DM, however, is the presence 
of an annual modulation in the event rate at direct detec- 
tion experiments [T] |2] . These experiments aim to detect 
the energy deposited by a nucleus after scattering with 
DM. The modulation signal arises due to the motion of 
the Earth relative to the galactic halo and peaks when 
the Earth is travelling fastest with respect to the halo. 
For most dark matter halo models, this occurs in late 
May or early June (see e.g. |3ll5). 

Two experiments have now detected an annual mod- 
ulation signal that has many features consistent with 
a signal arising from DM. The DAMA collaboration 
[SI [6] has observed an annual modulation over thirteen 
years; first with the DAMA/Nal setup and later with 
the DAMA/LIBRA upgrade. The large mass of the tar- 
get material and the long exposure time mean that the 
statistical nature of the modulation is beyond question. 
More recently, the CoGeNT collaboration 7J has pre- 
sented tentative evidence for an annual modulation after 
analysing fifteen months of data. 

In this paper, we assume both signals arise due to 
DM and study the consistency of the two signals as- 
suming the DM scatters elastically. In particular, we 
consider whether the phase and amplitude of the mod- 
ulations are consistent. Previous studies [SHU] have 
shown that DAMA and CoGeNT can be in agreement 
for DM of mass ^ 5-14 GeV that scatters elastically via 
a spin-independent interaction, albeit with tension from 
the CDMS J^, XENONIO il4j, XENONIOO ^ and 
SIMPLE [TB] experiments."'^ However, it is well known 
that interpreting the signals at direct detection experi- 
ments is sensitive to many uncertainties, particularly un- 
certainties in astrophysical parameters (see e.g. 



Therefore, we follow the approach of |27j, which specifies 
how to directly map experimental signals from one detec- 
tor to another, allowing a comparison without specifying 
any astrophysical parameters. Given the modulation ob- 
served at CoGeNT, we calculate the peak day and modu- 
lation amplitude expected at DAMA, which we compare 
with that observed at DAMA. The phase and amplitude 
will be the same if both modulations arise from DM. We 
finish by considering constraints from the low threshold 
analysis of CDMS II [131 [28] and XENONIO [l4|. 



II. DARK MATTER DETECTION 

We briefly review the formalism behind direct detec- 
tion experiments and define our notation. The differen- 
tial event rate for elastic spin-independent DM-nucleus 
scattering in counts per day per unit nucleus mass per 
unit exposure time per unit energy (cpd/kg/keV) as a 
function of recoil energy E'^; is 



dR 

dER 






F^{ER)e{ER)g{v^i^,t). (1) 



Here m^ is the DM mass, p^ the local DM density, /x„^ 
the DM-nucleon reduced mass and (T„ is the DM-neutron 
cross-section at zero momentum transfer in the elastic 
limit. F{Er) is the form factor which we take as 
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with q = ^/2mNER, R = J c^ + ^Tr^a^ — Ss^, 

c = 1.23^1/3 - 0.60 fm, s = 0.9 fm, a = 0.52 fm and 
niN is the mass of the target nucleus. We define Ct as 



CT = ti{UZ + UA-Z)f 



(3) 
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^ See I17fll9| for a critical discussion ot these experiments. 



where k, is the detector mass fraction of the target nu- 
cleus, A and Z are the nucleon and proton number re- 
spectively, while /„ and fp encode the coupling of the 
DM to neutrons and protons respectively. e(i?_R) is the 
efficiency of the detector and, in general, depends on Er. 



Finally, g{vmimt) encodes all the information about 
the DM velocity distribution: 
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where /locai is the local DM velocity distribution. Wmin is 
the minimum DM speed required for a nucleus to recoil 
with energy Er: 
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where fi is the reduced mass of the DM-nucleus system 
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III. A FORMALISM FREE FROM 
ASTROPHYSICS 

Following the approach of [27], we wish to map the 
differential event rate into Wmin space. As we will show, 
when the Umin space probed by two different experiments 
is the same, we can compare them directly without mak- 
ing any assumptions about any astrophysical parameters. 
This allows us to map the peak day and amplitude of the 
CoGeNT modulation onto the DAMA detector, which 
can be compared with the signal that DAMA observes. 
Below we briefly recap the relevant theory from j27) . 

Using Eqn. [5] we rewrite the event rate in Umin space: 
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where E^ = 2^'^Vmin/'m'N is now a function in Wmin space. 
t{Ep) and F{E^) are almost flat across the range of re- 
coil energies we consider at DAMA, CoGeNT and later 
CDMS II. Therefore, we remove them from the integral 
and evaluate them at the average value across the domain 
of integration. The error introduced by this approxima- 
tion is less than a few percent. 

All of the physics that determines the peak day 
of the modulation is encapsulated in the integral 
/ '"*'' vg{v,t)dv. By arranging to have the same integra- 
tion limits for DAMA and CoGeNT, the respective mod- 
ulations will have the same peak day if they both arise 
from DM scattering. Given an energy range at DAMA 
(D), we can invert Eqn. [s] to find the energy range at 
CoGeNT (C) which spans the same Wmin space. Doing so 
we find: 
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Given the modulation amplitude AR = (i?(tmax) ~ 
-R(^min))/2 observed by CoGeNT (in cpd/kg), and a spe- 
cific choice of my, and Ct, we can invert Eqn. [Tito solve 



FIG. 1: The peak day of the modulation signals measured by 
DAMA (red) and CoGeNT (blue) for different values of m^. 
At DAMA, we fit to the data in the 2-6 keVee energy range 
measured over thirteen annual cycles. At CoGeNT, we fit to 
the data in the energy range determined using Eqn. [8] It is 
clear that the GoGeNT modulation generally peaks earlier in 
the year than DAMA's modulation signal, although there is 
agreement at la for m^ ~ 7 GeV. The Standard Halo Model 
predicts the modulation peaks on Day 152 = 2nd June. 



for J"'""'' vg{v,t)dv. This allows us to calculate the am- 
plitude expected at DAMA (in cpd/kg): 
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We compare this with what is observed at DAMA. If 
both modulations arise from DM scattering, the expected 
and observed amplitude will be the same. Since no as- 
trophysical parameters enter into this formula, the two 
experiments can be compared free from astrophysical un- 
certainties. 



IV. COMPARING DAMA AND COGENT 

For the purposes of this analysis, we assume the DM 
only scatters off the sodium nuclei at DAMA. This is a 
very good approximation for the light DM that we con- 
sider. DAMA presents their time-dependent results in 
fixed energy ranges: 2-4 keVee, 2-5 keVee and 2-6 keVee. 
This can be converted to keV by dividing by the quench- 
ing factor for sodium, which we take as (jNa — 0.3. Vary- 
ing gNa does not lead to significant differences in our 
results. We do not include channelling, as theoretical 
calculations indicate this is small [29 . 

The CoGeNT collaboration has released the time- 
stamped data for independent analysis. We bin their 
data in the energy range that spans the same Wmin space 
as DAMA. Statistics at CoGeNT are limited, there- 
fore we concentrate on the 2-6 keVee bin at DAMA, 
as this corresponds to the largest energy range at Co- 
GeNT and will therefore contain the most events. To con- 
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FIG. 2: Best fit regions for DAMA and CoGeNT from a two parameter fit to the amplitude and peak day of the respective 
modulation signals. Day 152 = 2nd June. The blue region shows the amplitude and peak day expected at the DAMA detector 
based on what CoGeNT observe. The red contours show the regions obtained from fitting to the 2-6 keVee DAMA data in 
[5l[6]. In the left (right) panels m^ = 7 (12) GeV. In the upper (lower) panels fn/fp = 1 (—0.7). The purple solid line indicates 
the unmodulated rate at CoGeMT. Modulation amplitudes above this line are excluded. Note the lower panels have a different 
scale for the modulation amplitude. 



vert from keV to keVee at CoGeNT, we use i?[keVee] — 
0.199ii^[keV]^'^^ [30]. For each energy range, we divide 
the 458 days of data into 15 time bins with a width of 30 
days, and one final bin with a width of 8 days, correcting, 
when appropriate, for periods when the detector was not 
running. 

In Fig. [T] we show the peak day found at DAMA and 
CoGeNT. The regions between the red lines indicate the 
preferred regions at la (sohd), 90% (dashed) and 2a (dot- 
ted) after fitting to the residual events as a function of 
time measured by DAMA/Nal [5 and DAMA/LIBRA [5] 
over thirteen annual cycles. For each mass, we bin the 



data at CoGeNT in the relevant range determined from 
Eqn. [8] and plot the results in blue. For both DAMA 
and CoGeNT, we fix the period of the modulation to be 
one year and consider variations of the phase that sat- 
isfy x^ < Xmin + ^X^ ■ The variation in the peak day for 
different values of m^ is expected, since for each mass 
the data is binned in a different energy range. As well 
as purely statistical fluctuations in the observed count 
rate, N-body simulations have shown that variations in 
the peak day as a function of energy should be expected 
P2] and are also observed by DAMA (see Fig. 9 of [6 ). 
Generically, the CoGeNT modulation peaks earlier in the 



year than DAMA. It is only for m^ ^ 7 GeV that there 
is agreement at la. 

Fig. [2] shows the la (soUd), 90% (dashed) and 2a (dot- 
ted) contours in the peak day against modulation ampli- 
tude plane for m^ = 7 GeV (left panels) and m^ — 12 
GeV (right panels). In the upper panels, we choose 
fn/fp = 1 while in the lower panels fn/fp = —0.7. The 
choice fn/fp = —0.7 is phenomenologically motivated 
by the desire to suppress the event rate at xenon ex- 
periments [HI [5TH5i] . Given the modulation amplitude 
measured at CoGeNT, we show in blue the preferred re- 
gion calculated using Eqn. |9]for the expected amplitude 
at DAMA. The red region show the fit to the DAMA 
data measured over thirteen cycles. In our fits, we fix 
the period to be one year and subtract a constant rate 
at CoGeNT, determined from the value which minimises 
the x^- We proceed to find the values of the modulation 
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amplitude and peak day which satisfy x ^ xi 
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The energy ranges at CoGcNT for m^ = 7 GeV and 
m^ = 12 GeV, determined from Eqn. [8) are 0.70 - 2.38 
keVee and 0.87 - 2.96 keVee respectively. In these en- 
ergy ranges, the presence of an annual modulation is 
preferred over a constant event rate at 2. Oct and 2Aa 
respectively. The purple solid horizontal line indicates 
the constraint from the unmodulated rate measured by 
CoGeNT. We integrate the unmodulated rate in the en- 
ergy range, calculated using Eqn. |8j after subtracting 
the L-shell EC contribution and a constant background. 
We map this onto the DAMA detector using Eqn. [9J A 
modulation amplitude above this line predicts a modu- 
lated rate which is larger than the unmodulated rate. A 
modulation amplitude above this line is excluded since 
it would predict a negative event rate in the winter. Al- 
though the regions favoured by CoGeNT are large due 
to the low statistics, we can already see that m^ ~ 7 
GeV is preferred over 12 GeV and that the fit is worse 
for fn/fp — —0.7, where a larger modulation amplitude 
is generally predicted. 

For 7 GeV and fn/fp = 1 (upper left panel), we see 
that there is good agreement between the DAMA and 
CoGeNT modulation amplitude and peak day. However, 
the modulation amplitude is close to the unmodulated 
limit, so a large modulation fraction ('^ 70%) is required 
to be consistent. We will return to this issue in the next 
section. For fn/fp = —0.7 (lower left panel), the best 
fit modulation amplitude predicted by CoGeNT is much 
larger than that observed by DAMA. However, the 90% 
regions do overlap and a smaller modulation fraction (~ 
15%) is required to be consistent with the unmodulated 
rate. 



For 12 GeV and fn/fp = 1 (upper right panel), the 
CoGeNT unmodulated rate excludes all of the DAMA 
region. For fn/fp = —0.7 (lower right panel) it is only 
the DAMA and CoGeNT 2a regions that are in agree- 
ment. We thus conclude that a 12 GeV DM particle 
with fn/fp = 1 or —0.7 is disfavoured, before applying 
constraints from other experiments. 
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FIG. 3: The 90% lower confidence limit on the fractional mod- 
ulation required to be compatible with CDMS II. The blue 
(green) lines indicate constraints from the Ge (Si) analysis. 
The dotted (solid) line is for fn/fp = 1 (—0.7). For compari- 
son, the purple horizontal dashed line shows the modulation 
fraction from the Standard Halo Model (SHM). 



V. OTHER CONSTRAINTS 

CDMS II has analysed data over a similar energy range 
as CoGeNT in its low threshold Soudan Underground 
Laboratory germanium analysis |13| and Stanford Un- 
derground Facility silicon analysis |28j.^ Limits on an 
annual modulation signal have not been published but 
we can restrict the unmodulated rate, which we use to 
constrain the modulation fraction 



R{t, 
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For the germanium analysis, we use the 35 kg-days of 
data collected between October 2006 and September 2008 
from the T1Z5 detector, which has the best ionisation res- 
olution. For the silicon analysis, we use the 24.64 kg-days 
collected between December 2001 and June 2002. To set 
conservative limits, we assume that all of the observed 
events arise due to DM. In Fig. [3J the 90% lower con- 
fidence limits on the modulation fraction for fn/fp — 1 
(dotted) and fn/fp = —0.7 (solid) are shown. The blue 
(green) lines are from the germanium (silicon) analysis. 
We assume the modulation amplitude at DAMA is 0.0028 
cpd/kg/keV, which we see from Fig. [2J is on the lower 
edge of the CoGeNT 90% region for m^ = 7 GeV. For 
each value of the mass, we use Eqn. |8]to find the energy 
range at CDMS II which spans the same Wmin space as 
the 2-6 keVee range at DAMA. As the energy range varies 
for each mass, discrete jumps occur whenever a measured 



^ Silicon data has also been presented in |35l . We do not consider 
it here due to uncertainties in calibrating the energy scale for 
nuclear recoils near threshold 1361. 



event at CDMS II enters or leaves the energy range. For 
comparison, the purple dashed horizontal line shows the 
modulation fraction predicted from the Standard Halo 
Model (SHM). We see that the modulation fraction re- 
quired at DAMA and CoGeNT to be compatible with 
CDMS Il-germanium is larger for both choices of /„//p, 
while for CDMS Il-silicon, is larger for fn/fp — —0.7. 

The 2-6 keVee range at DAMA corresponds to 1.84- 
5.52 keV at a xenon-target experiment (for m-^—l GeV). 
We do not consider hmits from XENONIOO JS, as C^s 
has not been measured below 3 keV [37]. However, 
we can apply the limits from the S2 only analysis of 
XENONIO [11], which has a low energy threshold of 1.4 
keV. The 15 kg-days of data were collected between 23rd 
August and 14th September. Since this is approximately 
half way between the maximum and minimum of the Co- 
GeNT and DAMA modulation signals, we assume the 
measured rate is the same as the unmodulated rate. We 
use the parameterization for Qy given in [141 and assume 
all events arise due to DM. We do not apply the edge (in 
z) event rejection. For each mass, we again use Eqn. p^ 
to find the energy range that spans the same v^^^ space 
as the 2-6 keVee range at DAMA, and assume a modu- 
lation amplitude of 0.0028 cpd/kg/keV. For fn/fp = 1, 
we find the required modulation fraction is > 100% for 
all masses we consider. Hence, under the assumptions 
we have made, XENONIO excludes all parameter space. 
For fn/fp = —0.7 the constraints are severely weakened; 
a modulation fraction greater than ^ 2% is required over 
the whole mass range, which is easily achieveable. 



VI. SUMMARY AND CONCLUSIONS 

We have presented a comparison of the CoGeNT and 
DAMA modulation signals free from astrophysical uncer- 
tainties, having assumed that both modulation signals 
arise due to elastically scattering DM. 

We found that the peak day of the CoGeNT modula- 
tion is always earlier than the peak day of the DAMA 
modulation. The la confidence regions do overlap in 
some parameter space, but the best fit points typically 
differ by ^ 30 days. The SHM, which assumes the DM 
is distributed isotropically, predicts a peak day close to 
2nd June. If CoGeNT continue to measure a peak day at 
the lower edge of the DAMA confidence regions (around 
mid-May), there will be interesting consequences for DM 
galactic halo models. 

For DM that couples equally to protons and neutrons 
{fn/fp = 1), we found the measured modulation am- 
plitude and peak day at DAMA is consistent with that 
expected based on the CoGeNT results. However, the 
XENONIO S2 analysis excludes the whole mass range in 
question. Even if we were to ignore the XENONIO analy- 



sis, tension remains with the unmodulated CoGeNT rate 
and the constraint from CDMS Il-germanium, which re- 
quire large modulation fractions. For m^ = 12 GeV, the 
CoGeNT unmodulated rate excludes the DAMA modu- 
lation signal. Moreover, for all masses we consider, the 
modulation fraction needs to be larger than ^ 70% to 
be consistent with the low energy analysis of CDMS Il- 
germanium. The SHM predicts ~ 9%, which is signif- 
icantly smaller. Such large deviations from the SHM 
seem unrealistic. Therefore, based on the constraints 
from XENONIO, CDMS II and the CoGeNT unmodu- 
lated rate, we conclude that elastically scattering DM 
with fn/fp = 1 is unlikely to be the source of the DAMA 
and CoGeNT modulation signals. 

We also considered DM with isospin violating cou- 
plings fn/fp — —0.7. For this choice, the expected modu- 
lation amplitude at DAMA, calculated from the CoGeNT 
measurement, is generically higher than that observed at 
DAMA. For m^ = 12 GeV, we find that only the 2cr Co- 
GeNT and DAMA regions overlap. However, for m^ = 7 
GeV, amplitudes at the lower end of the CoGeNT 90% 
region are compatible. Furthermore, the CDMS II ger- 
manium and silicon constraints on the modulation frac- 
tion are milder, typically requiring modulation fractions 
^ 15%. This is still larger than that from the SHM, so 
if the DAMA and CoGeNT signals do arise from DM, 
the constraints from CDMS II indicate that there will be 
interesting consequences for the DM galactic halo model. 

It is clear that the evidence for a modulation in the 
CoGeNT data is still tentative and that much more 
data is required to definitively confirm a modulation. 
Fortunately, the CoGeNT-4 upgrade should provide 
much more data and significantly shrink the CoGeNT 
best fit regions. With more data, the approach presented 
here will serve as a useful complementary test on the 
consistency of DAMA and CoGeNT. In comparison to 
the usual method of displaying results in the (T„ - m^ 
plane, this approach has the advantage that the results 
do not depend on any astrophysical parameters. 

Note added: Following the submission of v.l of this 
work to the arXiv, related work appeared: [31] and ;39j . 
Our conclusions are in agreement. 
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